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We describe procedures to obtain the electronic structure of disordered systems using either tight 
binding hke models or quite directly from ab inito density functional band structure calculations. 
The band structure is calculated using super cells much larger than those containing a single minority 
component atom. We average over a large number of different super cell calculations containing 
different randomly positioned minority component atoms in the super cell as well as a varying 
total number of minority component atoms, weighted by the statistical probability. We develop an 
efficient and simple algorithm for unfolding of these bands, based on Bloch's theorem. The unfolded 
band-structure obtained in this way exhibits momentum and energy broadened structures replacing 
the gaps observed in often used single super cell calculations. Using the super cell averaged band- 
structure one can introduce a self-energy, resulting from the scattering of randomly positioned alloy 
components. The self-energy is causal, and shows strong energy and some momentum dependence. 
The self-energy shows rather non-trivial behavior and is in general non-zero at the Fermi-energy, 
resulting in an ill or undefined Fermi surface. The real-part of the self-energy at the Fermi-energy 
relates to an apparent violation of Luttinger's theorem. There is no simple relation between the 
apparent Fermi-surface volume and the electron count. Examples introducing these effects both for 
model and real binary alloy systems are presented. 

PACS numbers: 71.15.-m, 71.23.-k, 79.60.Ht, 71.20.Gj, 71.20.Be, 71.20.Eh 



The material specific properties of many systems de- 
pend crucially upon chemical substitution or doping, i.e. 
the intentionally introduction of impurities in an other- 
wise pure crystal structure. Substitutions can be cho- 
sen to introduce or change, the number of charge carri- 
ers, local magnetic moments, specific optical properties, 
or vortex pinning centers, just to mention a few com- 
monly known applications. These substituted atoms also 
introduce a destruction of translational symmetry and 
therefore scattering into the band structure as an im- 
portant deviation from the pure system. The effects of 
disorder resulting from random substitutions alter ma- 
terial properties in specific ways and the vast field of 
materials applications is strongly dependent on random 
substitutions.-^ -^^ 

Although we have powerful ah initio band structure 
methods to treat pure transitionally symmetric systems 
at least for systems in which electron correlation effects 
are less important these methods are generally not ap- 
plicable directly to randomly substituted systems or dis- 
ordered alloys. There are several useful approximation 
schemes available to treat randomly substituted systems 
or alloys, which can treat the disorder to a high level 
of accuracy. ^^^^ It is desirable to create a method able 
to describe disordered systems with the same level of un- 
derstanding and approximations as one has developed for 
periodic systems. Density Functional Theory (DFT)^, 
within the Local Density Approximation (LDA) is an ab- 
initio method proved to be useful for a large class of ma- 
terials. As long as the materials under consideration are 



not too correlated, the one electron wave functions in 
a Kohn and Sham^^ implementation can account for the 
material specific band-structure. First principle methods 
used for the description of crystals use the periodicity of 
the system in order to make calculations tractable. Two 
distinct problems arise when DFT is used to study the 
electronic structure of real materials: 1) how to account 
for disorder breaking the translational symmetry and 2) 
how to present the results in a meaningful way. 

Often, ab-initio calculations including disorder are 
done assuming periodic boundary conditions within a 
so-call super cell approximation, i.e. using a cell which 
is a multiple of the primitive unit cell. It is also typi- 
cal to choose the smallest cell size possible, containing 
one impurity atom, in order to reduce computational 
costs. Of course, the smallest cell with only one impurity 
or dopant, combined with periodic boundary conditions, 
can not account for a random impurity distribution in 
real materials. Such calculations basically result in the 
band structure of a fictitious new perfectly ordered ma- 
terial. Although in some cases such as the extreme dilute 
limit where the impurity impurity interference effects can 
be neglected this may be a reasonable approximation to 
the real problem it is important to develop theoretical 
methods based on the same approaches and approxima- 
tions used for pure systems to check for the importance 
of randomness. For a detailed discussion of development 
in ab-initio theory of alloys we refer to a recent review 
paper by Ruban and Abrikosov.^^ We too use the super 
cell approximation but allow for randomness in the impu- 
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FIG. 1: (color online) (a) Tight binding model for an infinite 
chain of dimers with intra-dimer hopping t and inter-dimer 
hopping parameter t^ (b) Band-structure of the model de- 
scribed in panel (a) with t' = t and the fatness of the bands 
calculated according to Eq. (O scaled by a factor of 0.2 for 
clarity. The doted lines are the dimer BZ. (c) same as panel 
(b) but now for (from left to right) t' = 0.75t, t' = 0.50t, 
t' = 0.25t, t' = O.OOt. 



rity distribution across the super cell. We then propose a 
simple procedure to present the unfolded calculated band 
structures in a way as seen by angle resolved photo emis- 
sion spectroscopy (ARPES) and demonstrate the effect 
of impurity scattering on the band structure and the self 
energy. 

We first present the method used to unfold the band- 
structure in the first Brillouin zone of the super cell to 
the first Brillouin zone of the original cell. Which is ba- 
sically based on Bloch's theorem and symmetrization of 
the wave- functions. The method can be readily applied 
to all first-principle codes able to calculate partial den- 
sities of states and local symmetry or orbital projected 
band structures. We then demonstrate the basic con- 
cepts of the random impurity calculations using a simple 
tight binding like model and on site impurity potentials. 
In the last section before the conclusion we describe the 
use of these concepts in the DFT calculation of a real 
alloy system composed of Li substituted with either H or 
Mg. 



I. UNFOLDING AND RANDOMNESS OF 
IMPURITIES 

In order to recover the original band structure from 
a super cell calculation one needs to unfold the bands 
back to the first Brillouin zone of the original cell. 
Several methods are available, which either look at 
the momentum dependent Hamiltonian,^^ or at the 
eigenfunctions.^^ In case of some additional distortion, 
i.e. the super cell is not a perfect repetition of the orig- 
inal primitive-cell unfolding can become a rather non- 
trivial task. A possible method to define unfolded bands 
is to assign a certain weight to each eigenstate so that 
the eigenstates with a weight of 1 represent the original 
band structure while the eigenstates with a weight of zero 
identify the folded bands. In this section we first demon- 
strate how this can be accomplished with a rather simple 
procedure, after the example we present a more thorough 
derivation. 

Consider a one-dimensional tight binding model with 
one s orbital per site and a nearest neighbor hopping inte- 
gral t. The solution to this problem is the wave function 
that is symmetric at the zone center and antisymmetric 
at the zone boundary: 

1 ^ 

E(k) = 2tcos{ka), 

where t < 0. 

As a super cell we take a doubling of the unit cell. This 
could model an infinite chain of dimers with intra-dimer 
hopping t and inter-dimer hopping parameter Assum- 
ing equal on-site energies and constant inter-atomic dis- 
tances of size a, the matrix representing the Hamiltonian 
is given by: 

/ t + t'e-^^>^-\ . 

The Hamiltonian can be understood by looking at Fig. 
([1]). There are two sites within the unit-cell, labeled A 
and B. (The unit cell, indicated by the black box has 
length 2a.) Starting from the left site in the unit-cell 
at the atom labeled A one can hop to the right with 
the strength t or to the left with strength f. Both pro- 
cesses represent a hopping from site A to site 5, and 
thus enter as off-diagonal elements between (j)A and (j)B- 
In the Hamiltonian the hopping is multiplied by a mo- 
mentum and directional dependent phase, whenever the 
hopping connects two species in a different unit-cell. For 
a hopping crossing the super cell boundary the additional 
phase is e~*^^ for a hopping to the right, or e+*^^ for a 
hopping to the left. With R = 2a, the size of the super 
cell. For t = t' the eigenvalues and the eigen-vectors are: 

m = \ll{4>A±e''"'4>B) (3) 
E{k) = ±2tcos{ka), 
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The main band has the energy of +2t (remember that 
t < 0) at the zone center and —2t at the zone boundary. 
The wave-function of the unfolded band within the super 
cell consists of the sum of (j)A and at the zone center 
and the difference at the zone boundary {k = n/a). The 
second solution, i.e. the folded replica, shows the oppo- 
site behavior in energy and its wave-function at the zone 
center (boundary) is made by the difference (sum) of (j)A 
and (j)B- Both bands can be seen in panel (b) of Fig. ([1]). 

In a true multi-component system one often uses a 
so-called 'fat' band representation in order to separate 
various orbital contributions into the band structure at 
given energy and k-point. There the weights (i.e. k de- 
pendent band thickness) are determined by the squares 
of the eigen-vector components: 
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FIG. 2: (color online) (a) LDA Li band- structure calculated 
for the primitive unit cell, (b) Li band-structure calculated 
for a 2 X 2 X 2 super cell. In both panels colored fat bands are 
used to show the 2p and 2s orbital character within the first 
Brillouin zone of the primitive unit cell. Unfolding is done 
with the use of Eq. (O, traced over the 2p and 2s orbital 
character. The zero of energy is at the Fermi energy. 



with b the band index, r and some orbital or site in- 
dex, C the crystal momentum {k) dependent eigen-vector 
component, and i a sum over equivalent sites within the 
unit-cell. It is clear that weights calculated using the 
eigen-vectors and weights as defined above results in a 
constant weight throughout the whole Brillouin zone for 
both bands. 

A simple change in the definition of the weighting- 
function allows one to include information about the rel- 
ative phase in the wave- functions between different sites. 
If one defines the weight as: 




(5) 

one basically changed the sum of the norm-squared of 
the eigenstate pre-factors into the norm of the Fourier 
transform of the eigenstate pre-factors. Defining the 
positions of the atoms such that = and tb = 
a and applying Eq. ([5j) to the wave-functions as 
found in the dimer model immediately shows that the 
band belonging to t he o riginal, unfolded system, with 
the wave- function y^^l/2{(j)A + g^^^^b) gets the weight 
{l/2){^/l/2 ^ {e-'^''e'^'')^yT/2)^ = 1. The folded band, 
with the wave- function y^l/2{(j)A — e^^'^cjys) obtains the 
weight, {l/2){^lj2- {e-'^''e'^'')^lj2Y = 0. 

As a proof of principle, we show in Fig. ([2]) the band 
structure of Li metal calculated with TB-LMTO density 
functional code^^ using a 2 x 2 x 2 super cell. The intensity 
is calculated according to Eq.(|5]) in both super cell and 
primitive cell calculations. Li is a sp metal with a body 
centered crystal structure and a lattice constant of 3.51 
A. The primitive-cell contains one Li atom, leading to 4 
bands per /g-point of different orbital character. At F, 
for example, the Li 2s derived band is at ^-4eV. The 
three fold degenerate 2p bands are at about 14 eV above 
the Fermi energy. In other regions these bands cross and 
intermix. The variable r labels the orbital character: 
r G {2s^2px^2py^2pz}. The index i sums over the 16 



equivalent sites within the super cell. The band-index 
h labels the 4 x 16 different eigenstates found at each k 
point within the super cell. To highlight the complete 
band structure, we use the total intensity, traced over 
the orbital index r: ll^^ = I^'^. 

There is a complication with the definition for unfold- 
ing presented so-far. In principle one is free to choose the 
phase of each basis function separately. Multiplying the 
phase of the basis function in the dimer labeled as A in 
Fig- ([U) should not change any physical results. How- 
ever it does change the unfolded weight as calculated 
with Eq. (j5j). The solution is simple, one should keep 
track of the phases of the basis functions and introduce a 
complex- conjugated of this phase in the definition of Eq. 
(|5j). Below we will derive the equation for the weight of 
an unfolded band-structure more rigorously. The here 
sketched phase-problem will then come out naturally as 
additional form-factors in the sum of Eq. ([5]). 



A. Derivation of general unfolded weight of a band 

The weight of a band can be calculated in similar fash- 
ion as one would calculate the momentum (inter-, and 
intra-Brillouin-zone) depended photo-emission spectra. 
Neglecting polarization effects, assuming a high photon 
energy such that the final states can be approximated by 
plane waves, and neglecting the fact that the action of 
the electric field of the light on the electrons is given by 
a momentum operator one can write the spectral weight 
(/) of a one particle eigenstate ip as: 



4,k cx (e*^-^|V^5,fc(r)) 



(6) 



whereby k labels the momentum, k the corresponding 
reduced or crystal momentum, and b the band-index. (A 
momentum dependent mixture of 25, 2px, 2py^ or 2pz 
orbital character in the case of Li as shown in Fig. (|2]).) 

For a super cell calculation of an undistorted structure 
it is straight-forward to show that the spectral weight 
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function as defined in Eq. (j6j) will be zero for bands that 
do not belong to the unfolded bands of the primitive unit 
cell. Take the primitive cell to be defined by the matrix 
A, consisting of three vectors ai, a2, and as, which span 
the parallel piped defining the primitive unit cell. Using 
Bloch's theorem for a periodic system we know: 



(7) 



with, 



Ub,k{r + n • A) = ut,k{r), Vn G Z^^^. (8) 



Here n is an integer vector {n - A = niai + n2a2 + tisSls) 
and is a periodic function with the period of the 

primitive unit cell, ?/^5^fc(r) is periodic up-to a phase: 



^/;fc,fe(r + n-A) = e''=-(»-^)VM(r)- 



(9) 



For a calculation done in the primitive unit cell one can 
use the block diagonal in k characteristic of the Hamil- 
tonian. Solving H for each k vector separately yields the 
band-structure. Within a super-cell calculation one does 
not obtain the momentum of a wave-function directly. 
If one takes a super-cell without perturbing the Hamil- 
tonian the eigenstates are non-the-less given by Bloch's 
theorem as given in Eq. ([7]). The Fourier transform, 
of such a function, which defines the spectral weight as 
defined in Eq. ([6]) is given as: 



(10) 



= I E / e-''' (--+"^ -^)tA6,fc(r + n,- ■ A) dv 
I ? In, 



= IE' 



Ak-k)-(n,-A) 



[ e^(^-'')--M6,fc(r)ar 



with k the crystal momentum of the wave-function in 
the primitive unit cell. The reduced momentum k^ is an 
eigenstate property according to Bloch's theorem, but 
undetermined from the super-cell calculation. The re- 
duced momentum index within the super-cell has been 
suppressed as it is up to a super-cell reciprocal lattice 
vector implicitly labeled by k. Q labels an integral over 
all space and ftp an integral over the original primitive 
unit cell, and the sum over all possible independent shifts 
by an integer lattice vectors {n ■ A) of the primitive unit 
cell. The first term of the unfolded spectral weight in the 
primitive unit cell (/6,k) contains a sum over all unit-cells. 
This term leads to the selection rule: 



k-k=nG, 



(11) 



with n an integer vector and G a matrix representing 
the reciprocal lattice: = 2ttA~'^. This term is re- 
sponsible for the fact that only those bands that belong 



to the principle unit cell will gain weight, i.e. it unfolds 
the folded band-structure. The second term of the un- 
folded spectral weight is a material and Brillouin zone 
dependent form factor. In photo-emission experiments 
it is this term (together with the here neglected polar- 
ization dependence) that is responsible for the Brillouin 
zone dependent spectral weight. 

The unfolding of bands as described in Eq. (j6j) can be 
simplified further by approximating the Brillouin zone 
dependent form factor. We will first briefiy describe 
model calculations and then extend it to a formalism 
that can readily, without much effort, be implemented 
in most ah initio density functional programs. 



1. Unfolding of hands in tight-hinding model calculations 

For model calculations one often expresses the band- 
structure in terms of some form of tight-binding repre- 
sentation. If we take the variable r as the index for the 
local orbital and or spin character of the tight-binding 
basis set (in the case of Li r is 2s, 2px, 2py, or 2pz) then 
one can define the spectral weight matrix for a given one- 
particle eigenstate {ipb^ki'^)) as: 



(V^6,fc(r)|aj^ ,a, |V^5,fc(r)), 



(12) 



whereby the operator ^ annihilates an electron at mo- 
mentum k and of orbital (spin) character r. The mo- 
mentum dependent annihilation operator can be written 
as a sum of position dependent operators as: 



1 



E 



(13) 



which describes the relation between the Fourier trans- 
form of a one-particle wave-function and the expectation 
value of the momentum dependent occupation number 
operator. If the orbital character of a certain band is of 
no importance one can write the total unfolded spectral 
weight of a certain band as the trace over the spectral 
weight matrix: 



^6,k 



(14) 



Unfolding of hands for ah initio density functional 
calculations. 



In the case of an ah initio density-functional theory cal- 
culation one could naturally first create a tight binding 
model and then use Eq. ([12]) in order to calculate the un- 
folded spectral weight. Such a scenario has recently been 
proposed by Wei Ku et al}'^ and successfully used on sev- 
eral systems ji^ii^ The unfolded band-structures found in 
this way are expected to be rather similar to the un- 
folded band-structure one obtains by defining the weight 
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of a band via a Fourier transform (see Eq. (j6j)). For 
large super cells with several impurity atoms, defining 
a tight-binding Hamiltonian might become troublesome. 
An alternative method that allows one to directly obtain 
the unfolded band-structure from an ab initio calculation 
therefore is useful. 

One can simplify the Brillouin zone dependent form 
factor of the unfolded spectral weight as it appears in Eq. 
([To]) by looking at the one-particle wave-function charac- 
ter within a certain (muffin-tin) sphere defined around 
each atom. These projections are implemented in most 
ab initio DFT codes, where these projections are used 
to calculate the partial character of a given eigenstate 
needed to plot partial density of states, or bands of a 
certain orbital character. The projected wave- function 
'06,fc(r) can be written as: 

V^6,fc(r) = -^'^Cli^i,(l)r,i{r - Yi), (15) 



The integral: 

/ e-^^->,,,(r)ar, (18) 

defines a Brillouin zone dependent form factor, and al- 
though this is a physical quantity one can simplify the 
unfolding of bands substantially by neglecting this term. 
If the orbital basis set is defined such that the main outer 
part of 0r,i(i*) is real and positive for all r then a power- 
ful approximation is to set the integral in Eq. (p!8|) equal 
to 1. In this case the unfolded weight becomes: 

^^k' = ]^(Ee"''^^"<a)*(Ee"'''''''C'^.,fc)- (19) 

i i 

In the case that the super cell is undistorted, i.e. a pure 
repetition of the system one can understand that this 
definition of the spectral weight indeed unfolds the band- 
structure with the argumentation based on Bloch's theo- 
rem. This can be done by following the argumentation as 
given in the beginning of this section from Eq. ([7]) to Eq. 
([To]). When the super cell is not a perfect repetition of 
primitive basis cells then the weight of the folded bands 
starts to deviate from zero. The weight in the folded 
bands is related to the mixing of the folded and original 
bands in the band-structure and therefore a measure of 
the distortion. 



with b a band index of the one-particle eigenstates, r 
an orbital and spin index, the sum over all sites i that 
contain a muffin tin sphere onto which one projects, 
the position of the sphere labeled with the index z, 
the functions projected onto (with orbital character r 
and centered at position r^), and ■ ^ the projected pre- 
factor. Similarly as done for model calculations it makes 
sense to define an orbital (r) specific unfolded weight, 
starting from: 



r,,u{r) = E Cl.,k4>r,i{r - ri), (16) 



Unfolding such a projected one-particle wave- function 
with the use of Eq. ([6]) yields: 



(17) 



One should be slightly careful how the form-factor as 
presented in Eq. ([T8|) is neglected. The phase of this in- 
tegral depends on the phase of the basis function ^^-^^(r), 
which can be chosen freely. In many calculations the pro- 
jections are done onto atomic like wave- functions. These 
functions are real, which fixes most of the phase, but 
still can either be positive or negative. The radial part 
of atomic like wave- functions has nodes. The important 
outer part of the radial wave- function is either positive or 
negative, depending on the number of nodes. The num- 
ber of nodes for an atomic like wave-function is given as 
n — l the phase of the basis function therefore is (—1)^^"^^. 
Naturally different implementations might use different 
phases for the basis functions. (A complex phase of 
g-*k ri example be a useful alternative choice.) 

Using our phase definition, including the factor (— 1)(^~0 
for the basis functions, the Brillouin zone dependent form 
factor will be approximated as: 

/ e-'^-^'^rAr) dr ^ (-1)(<-^^), (20) 

whereby nj is the principle quantum number of the or- 
bital with index r and position and IJ is the angular 
momentum of the orbital with index r and position r^. 



4?k' = (^e-'^-<,(r)ar)*(^e-'^'>,-,,(r)ar) 

= WS^L (r - ri) ar) * ( ^ e-^^-C,- (r - r,) ^r) 

= ]^ ( E e--^--^,, , (r) 5r) * ( ^ e-"^--C,-,,,, ^ e""^' (r) dr) . 

I 
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Distribution function for p=0.05 




FIG. 3: (color online) Statistical fraction of clusters with Ns 
substituted atoms for a cluster size of 4 (left top) to 128 (right 
bottom) Nt total atoms and an average substitution of 5% 
(y9 = 0.05). 



The unfolded weight then becomes: 

= ]^(E^"^''"(-i)'"-'"'''^<^)* (21) 

i 

x(E^"''"(-l)^"'"''^^M.^)- 



The partial character • ^ is the projected weight of the 
eigenstate with band index 6, orbital and spin character 
r, at site i in the super cell for the crystal momentum 
k. The unfolded weight of this band is calculated by 
taking a Fourier transform summed over (pseudo) equiv- 
alent sites {i) within the super-cell. This can either be 
the same atom at positions in the super cell slightly devi- 
ating from a perfect repetition of the primitive unit cell, 
or different atoms in the super cell in the case of sub- 
stitution. The norm-squared of this Fourier transform 
defines the unfolded weight of a certain band and orbital 
character. 



cluster is given by the binomial distribution: 



P(iV,,iV,;p)=p^«(l-p)W-^^) 



7V,!(7V,-7V,)!' 



(22) 



In order to illustrate the impurity distribution further we 
show in Fig. ([3]) a plot of this distribution for a doping 
of 5% {p = 0.05) and a cluster size of 4, 8, 16, 32, 64, and 
128 Nt total atoms. One should realize that even for a 
reasonable cluster size there still is a finite chance to have 
no substituted atoms within the cluster. For each total 
number of substituted atoms in the cluster there are still 
many different ways to arrange these impurities. It is 
important to take this into account, since the scattering 
of nearest neighbor impurities is in general very differ- 
ent from the scattering of isolated impurities. In general 
one should determine the band structure in a super cell 
method with as large a cell as possible in real space and 
average over a varying number of impurities and a vary- 
ing distribution of these impurities. Each of these calcu- 
lations should be weighted with the probability to find 
this particular kind of configuration, including both the 
chance to have this specific density within the cluster and 
this specific order within the cluster. In practice it is not 
needed to calculate all possible impurity configurations 
(including variation of density and distribution) within 
the cluster. A random subset will already give a reason- 
able approximation to the real solution, provided that 
the cell is large enough and a large number of random 
configurations are taken. 

We define a configuration to be a specific realization of 
a super cell with a random number of substituted atoms 
at random positions. The average density of substituted 
sites is given by p. 



II. RANDOM ONSITE ENERGY IMPURITY 
MODEL 



B. configuration average for the description of 
randomness 

There is a crucial difference between a calculation for 
an ordered structure and a system with randomly substi- 
tuted atoms. In a system with ordered impurities gaps 
will open in the band-structure, whereas for a random 
system this may not be the case and gaps can become 
smeared out structures, as will be shown by several ex- 
amples later. In order to simulate the situation of an 
infinite crystal with a random impurity distribution one 
would like to create very large super cells. In order to do 
these calculations with the use of tractable cluster sizes 
it is important to average over different random distribu- 
tions of impurities with equal impurity density and over 
different impurity densities within the cluster. Lets de- 
fine the chance to have a substituted atom at site i to be 
equal to p. For a cluster of Nt total sites the probabil- 
ity to have a number of Ns substituted atoms within the 



Here we exemplify several aspects of random impurity 
calculations, using a simple model system. The model 
we use is well studied which allows us to compare the su- 
per cell calculations to theoretical approximations used 
elsewhere.^-"— We use a single band model with con- 
stant nearest neighbor hopping (t) and a onsite impurity 
potential V at a fraction p of all atoms randomly placed 
within the solid. The Hamiltonian is given as: 

i7=^taJa^. + ^F4a,, (23) 

ihj) * 

with the sum over all sites i and nearest neighbor sites 
j. t is the hopping-strength, resulting in a band- width 
of W = 2Dt. D is the dimension of the system. Vi is a 
potential with the value F at a fraction p of randomly 
placed sites and zero at other sites, p is the density of 
substituted atoms. 

As the Hamiltonian only includes one-particle interac- 
tions it is natural to Fourier transform it to momentum 



7 



(k) space: 



H 



k 



-k^k^k 



En 



k,k''^k'^k'' 



(24) 



k,k' 



with Vk,k' = ;^ ^■'^^Fi and ri the position of 

atom i. Within this representation the hopping part be- 
comes diagonal in k. The random potential allows for 
scattering between different k- vectors. In the limit of 
large cluster sizes or for a potential averaged over many 
different cluster configurations one finds that Vk,k = Vp 
and that |Vk,k'| for k 7^ k' scales as ^/T/N. The aver- 
aged Hamiltonian (i^), is thus defined as = ^^^(^k + 
Vp)al^a^. One might be tempted to state that for in- 
finitely large clusters or for a cluster of finite size, aver- 
aged over a large number of configurations the Hamilto- 
nian becomes diagonal in momentum k and thus momen- 
tum is always a good quantum number. Care has to be 
taken as there is indeed only an infinitesimal small cou- 
pling to different momenta for infinitely large clusters, 
but there are an infinite number of different momenta 
possible to which one can scatter thus leading to a finite 
perturbation. 

The one particle Green function is defined as: 



Gk,k'(^) = lim { 

r^o+ 



1 



4|0). 



(25) 



The spectral function whose poles define the band- 
structure, is given as A(k, cj) = — ;^Im[Gk,k(^)]- The 
off diagonal elements of Gk,k'(^) in k, i.e. those ele- 
ments for which k ^ k' scale as a function of cluster 
size, in the same way as the off diagonal elements of 
the Hamiltonian, namely as y^l/N. For large clusters 
or for an average over many different random configu- 
rations the Green function becomes diagonal in momen- 
tum Gk(^) = Gk,k'(^)- Calculating the Green function, 
using the average Hamiltonian (which is diagonal in mo- 
mentum) is however different from using the full Hamilto- 
nian, even if one averages the Green function afterward. 
In the following whenever we write G\^{uj) we will assume 
that this is the Green function averaged over several con- 
figurations (Gk(c^) = Gk,k'(^))- 



A. Moment analyzes and series expansion 

The random onsite energy impurity model has been 
extensively studied with the use of series expansions. ^^^^^ 
The n-th moment (/i^^^) of the momentum integrated 
spectral function (density of states) is given as: 



(n) 



TV 



Tr[i^^] 



(26) 



with N the total number of sites in the cluster. The 
Hamiltonian is given as H'^^^ + i^(^) with = 

Xlk ^k<^k<^k ^^^^ ~ Xli ^i^i^i ' Analytical expres- 
sions can immediately be written down, if one neglects 



contributions to the moment which arise from products 
of i^^o) and H^^\ The moments of H^^^ are given as: 



(27) 



which are the moments of the original, unperturbed sys- 
tem. The non-zero lowest order moments are /ig^^ = 1, 
6t^D{2D - 1). The moments of 



,(2) 



/iQ ' = 2Dt\ and /ig 
are: 



(4) 



(n) 
Ml 



(28) 



Contributions due to products of H^^^ and H^^^ are im- 
portant for the moments with n > 3. 

The momentum k dependent Green function can be 
expressed in terms of a series in H^^^ with the use of the 
Dyson equations. For a Hamiltonian only containing one 
electron terms one can write: 

Gk.k'(c^) = ('^) + E (c.)<)„Gk".k'('^)- (29) 

k" 

For the random configuration averaged Green function, 
Gk,k'(^) = Gk(cj) one finds to first order in the onsite 
impurity potential V: 



Gk(w) = Gf (w) + Gf (w)ypGk(w). 



(30) 



It is often useful to define the self energy for the averaged 
Green function as: 



E(k,cj) 



1 



(31) 



If one compares the definition of the self energy with 
Eq. ([29|) one finds a striking similarity. Multiplying 

Eq. (jSU) with G^^\uj) and Gk(a;), one reproduces the 
form of Eq. (|29|) . but now for the configuration aver- 
aged, single momentum dependent Green function, in- 
stead of the original Green matrix. Using Eq. ([3Q|) one 
finds that to first order in V the self energy is given as 
E = Vp. To linear order in the perturbing potential one 
can thus replace the random impurity system by a pe- 
riodic system, which has an averaged potential on each 
site. This approximation is known as the virtual crystal 
approximation. It can correctly describe the averaged 
band-structure and total density of states to first order in 
V. It can not describe changes to the one-particle wave- 
functions important for the band-character, the spatially 
varying charge density, the broadening of bands, or the 
possible appearance of impurity bound states, which arise 
in quadratic and higher order of the perturbing potential. 



B. Cluster size and configuration averaging 

In order to fully describe the effects of random impuri- 
ties it is necessary to go beyond a series expansion in the 
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FIG. 4: (color online) Unfolded band-structure, density of 
states and Fermi- surface of a 2D tight binding model as pre- 
sented in Eq. (|23)) . for various cluster sizes. The left col- 
umn shows calculations for one random impurity distribution 
with p — 0.25. The right column shows calculations for an 
average of 1000 different random configurations, which only 
on average have an impurity density of 25%, and vary both 
in impurity density as well as impurity distribution within a 
cluster. In blue (red) the character of the host (substituted) 
atoms is shown. The potential of the substituted atoms is at 

V — — 0.3VK, with W the band-width (the potential of the 
host atoms is at 0). The Fermi-surface plots assume a filling 
of on average 1 electron per atom. 

perturbing potential. For such cases doing exact diago- 
nalization of large clusters is a powerful method, espe- 
cially with the current computational resources. Within 
this section we compare several calculations done for dif- 
ferent cluster sizes. We furthermore show that the clus- 
ter sizes can be considerably reduced if one averages over 
several random configurations, i.e. both varying the im- 
purity density as well as the impurity distribution within 
the cluster. 

In the left most column of Fig. (j4]) we show calcula- 
tions for the band-structure of a single random configu- 
ration, with fixed concentration. The different rows show 
the calculation for a different total cluster size. The cal- 
culations are done for a 25% (p = 0.25) substitution and 

V = — 0.30iy. The smallest cluster possible which allows 
for a substitution of 25%, namely with a total of = 4 
sites and one impurity atom, shows the opening of large 
gaps in the band-structure. This is to be expected as for 
a 2 X 2 cluster with 1/4 substitution one has a perfectly 
ordered super-structure, which is very different from a 
random system. The larger the cluster becomes the 
smoother the band-structure becomes. The gaps become 
smaller and smaller and more and more distributed to 



random places in momentum space. One however needs 
very large clusters to find convergence if one does not 
average over different random configurations. Even for a 
single calculation of a 32 x 32 super cell there is still some 
gaping visible in the band-structure (left column bottom 
row of Fig. dl])). This becomes much better if one aver- 
ages over different random configurations (i.e. averaging 
over both a varying number of substituted atoms within 
the cluster and averaging over the possible different ways 
to substituted these atoms). The forth column shows the 
same band-structure calculations, but now averaged over 
1000 random configurations. As one can see these look 
much more smooth. The gaps disappear when an average 
over different configurations is taken. Even for a 4 x 4 
super cell calculation the averaged band-structure looks 
quite reasonable, whereas the calculation of a single con- 
figuration is far from converged. What happens is that 
each different random configuration has gaps at different 
energies and k vectors, leading to an averaged broaden- 
ing of the bands and not to a gaping of the bands. A 
similar effect can be seen for the density of states. The 
calculations for a single configuration are rather noisy; 
the average over 1000 configurations however shows a 
nice continues behavior. 

Averaging is crucial for the calculation of the Fermi- 
surface. The Fermi-surface is calculated by assuming a 
filling of 1 electron per atom in the unit cell after av- 
eraging all calculations done for different configurations. 
A single configuration shows only some points on the 
Fermi-surface, whereas an average over several configura- 
tions shows a very nice, though broadened Fermi-surface. 
Note that in the Fermi-surface plots for small super cell 
calculations one can see the periodicity of the super cell 
very clearly as features with the periodicity of the super 
cell Brillouin zone. The use of small super cell calcula- 
tions can lead to spurious shadow bands. 



C. Comparison to other approximations 

Next we compare the direct method of calculating the 
band-structure, Fermi-surface and density of states as 
presented here with other methods like the coherent po- 
tential approximation and its non-local version.—*^ in 
Fig. ([5]) we present calculations as a function of the impu- 
rity potential for a 2D 50% {p = 0.5) substituted system. 
The model and parameter range used are the same as pre- 
sented in Fig. (5) of the paper by Batt and Rowlands^^, 
which allows for a direct comparison between the dif- 
ferent approximations. For a potential of V = —0.3W 
one finds a clear Fermi-surface at a filling of one electron 
per site. The Fermi-surface is rather broad. These find- 
ings are basically equivalent to the results found with 
the coherent potential approximation and for the non- 
local coherent potential approximation calculated using 
a 2 X 2 cluster. Increasing the impurity potential further 
will eventually lead to two separate bands separated by a 
gap. In principle one could state that the system becomes 
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of 100 random configurations, p=0.50, 2D 




function {G^\uj)) and the self energy as: 

1 



r M 



FIG. 5: (color online) Unfolded band-structure, density of 
states and Fermi- surface of a 2D model calculation as defined 
in Eq. ([23l), for a 16 x 16 cluster, with p = 0.50 and V ^ 
— 0.30VK to V = — 0.55VK. The Fermi-surface plots assume a 
filling of on average 1 electron per atom. 



insulating at half filling. However, due to the large broad- 
ening of the bands, one still finds for a very large potential 
range a finite intensity at the Fermi-energy. It is in this 
range that one finds differences between the coherent po- 
tential approximation and their cluster extensions^^^^^. 
For V = 0.55iy one finds a new Fermi-surface topology. 
Around F one finds an electron pocket which is related to 
the host band. Around M one finds a hole pocket derived 
from the impurity band. It is note-worthy that the co- 
herent potential approximation is not able to reproduce 
this effect. ^^^^^ The non-local coherent potential approx- 
imation, for large cluster sizes should naturally exactly 
reproduce the current calculations, under the condition 
that the cluster size taken into account is large enough. 
For a cluster size of 2 x 2 one can see that some features 
of the gaping are indeed correctly reproduced. The ex- 
act potential at which the Fermi-surface topology change 
takes place, as well as some details of the Fermi-surface 
are however different. For accurate calculations it is im- 
portant to take relatively large cluster sizes, within the 
non-local coherent potential approximation. The effect 
of the bath in the non-local coherent potential approxi- 
mation would be relatively minor for such large clusters 
sizes. It is therefore relatively straight forward to ex- 
tend the non-local coherent potential approximation to 
the current cluster calculations. For large cluster sizes 
both methods are equivalent. 



D. Self energy 

A powerful concept for the understanding of ran- 
dom impurity calculations is the self energy i;(k, cj) = 
— — . The self energy is a measure of how- 

much the band-structure of the system with impurities 
deviates from the system without substituted atoms. The 
full Green function is given in terms of the original Green 



S(k,w) 



(32) 



The real part of S defines an energy shift of Gk(c^), with 
respect to g\^\uj). The imaginary part of E is related to 
an energy broadening of the Green-function. The Self- 
energy is casual, meaning that in our representation the 
imaginary part is strictly negative, and goes to zero as 
l/o;^ for large absolute values of u. The real-part is re- 
lated to the imaginary part by Kramers-Kronig. Note 
that for a system of impurities there are no specific rela- 
tions for the self-energy at the Fermi energy. As seen in 
Fig. dH and (|5|) there is a substantial broadening of the 
Fermi surface (non-zero imaginary part of the self energy) 
as well as a possible change of Fermi surface topology (en- 
ergy and momentum dependent real and imaginary part 
of the self energy). 

It is generally excepted that a system with random 
impurities has a broadened band-structure compared to 
a pure system. The imaginary part of the self energy 
must be integrable and goes to zero for large and small 
energies. By virtue of the Kramers-Kronig relations this 
means that one then also must have a downward energy- 
shift for energies lower then the main broadening and an 
upward energy shift for energies higher then the energy 
at which the broadening takes place. The broadening of 
bands will thus go hand in hand with shifts in the band- 
structure. This will become more clear with the use of 
the example given in Fig. (j6j). 

In Fig. (j6j) we show the band-structure, density of 
states, Fermi-surface (for a filling of one electron per 
atom) and real and imaginary part of the self energy 
at r as well as their momentum dependence for a 31^ 
cluster as a function of the impurity potential at p = 0.3. 
The real part of the self energy at F (sixth column) has 
a constant part, independent of the energy, equal to Vp, 
which one would expect in linear order inV. In the forth 
column one can see the imaginary part of the self, energy, 
which peaks around V. This leads by Kramers-Kronig to 
a real-part of the self-energy, which is related to a shift 
of the band-structure. 

The shift of the band-structure due to the real-part of 
the self energy can be seen very clearly in the Fermi- 
surface plots. For a small potential of the order of 
V = — O.IOVF or smaller the band-structure, density of 
states and Fermi-surface are, aside from some broaden- 
ing, not changed much compared to the original host elec- 
tronic structure. The Fermi-surface of the host is shown 
in the third column, by the small green pockets around 
M. For V = — O.lOiy the Fermi-surface of the system 
with impurities overlaps with the original Fermi-surface. 
For a larger impurity potential the band-broadening be- 
comes more substantial, and a corresponding shift of the 
band-structure can be seen. For the Fermi-surface this 
results in apparently larger pockets around M. The elec- 
tron count does not change between the different cal- 



10 



Potential dependence, 10 x 10 x 10 cluster, average of 100 random configurations, P=0.30, 3D 




Im[£] (-Vp) Re[I] (-Vp) 

FIG. 6: (color online) From left to right; the unfolded band-structure, the density of states, Fermi- surface, the imaginary part 
of the self energy H at k = (F), the momentum dependence of the imaginary part of the self energy, the real part of the 
self energy E at k = (F), and the momentum dependence of the real part of the self energy. The calculations are done for 
the model as defined in Eq. (|23p . for a three dimensional 10 x 10 x 10 cluster, with p — 0.30 and the potential varies from 
V = —0.10 to y = —1.00 for the different rows shown. The Fermi-surface plots assume a filling of on average 1 electron per 
atom. The additional green small lines in the Fermi-surface plot show the fermi-surface for the unperturbed system with a 
filling of 1 electron per atom (small around M) and a filling of 0.7 electron per atom (larger around M). 



culations shown in the different rows in Fig. ([6]), The 
Fermi-energy is defined such that there are 1000 elec- 
trons in the cluster containing 10 x 10 x 10 atoms. It 
appears that Luttinger's theorem is violated. Calculat- 
ing the number of electrons in the system by looking at 
the Fermi-surface as defined as the maximum intensity 
of 74(k, cj) at the Fermi energy gives an incorrect num- 
ber of electrons. One can understand these findings by 
realizing that the calculations are not done for a periodic 
system, but averaged over many different randomly or- 
dered periodic systems, within a super cell. The Green 
function is obtained by unfolded these super cell calcula- 
tions back to the original Brillouin zone. The conditions 
for Luttinger's theorem to be applicable are not fulfilled. 

The change in Fermi-surface can be understood in a 
different way. For large negative impurity potentials the 
substituted atoms will build an impurity band, at a lower 
energy then the original host band. This impurity band, 
will be doubly occupied, i.e. full. This doubly occupa- 
tion of the impurity band removes p electrons from the 



Host band, leading to an apparent change in electron oc- 
cupation. In order to test this simple picture we added 
the Fermi-surface one would expect for an unperturbed 
system with p = 0.30 holes, i.e. on average only 0.7 
electrons per unit cell. These are the large Green Fermi- 
surface contours shown in the third column of Fig. (|6j). 
Indeed it seems that for large potentials the apparent 
hole doping can be understood in terms of doubly occu- 
pied bound states around the substituted atoms, which 
do not contribute to the Luthinger's count. There is no 
clear border between the two regimes where for small po- 
tentials the virtual crystal approximation is reasonably 
fulfilled to a regime for large potentials where the system 
can be understood in terms of an impurity band. There 
is always some deviation from the virtual crystal approx- 
imation, however only to quadratic order in the impurity 
potential. 
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FIG. 7: (color online) Calculations for: Left; a super cell of 
([100], [010], [002]). Middle; a super cell of ([110], [110], [001]). 
Right; a super cell of ([110], [101], [Oil]). Bottom row shows 
the super cell and the total self-consistent LDA energy with 
respect to the most stable arrangement. Middle row shows the 
folded and un-folded LDA band-structure of pure Li within 
the super cell. The bands with weight zero (should disappear 
in the unfolding scheme) are plotted in orange, the unfolded 
bands (with weight 1) in Black. Top row shows the unfolded 
bands for the Lio.75Mgo.25 ordered alloy within the 3 different 
super cell configurations. 



III. AB INITIO RESULTS 



There is an important difference between model calcu- 
lations and self-consistent first principle studies. Within 
a model the potential of the host and impurity are fixed, 
irrespectively of their electron count. Within a self- 
consistent calculation the potential depends on the local 
electron count. This later effect will lead to screening, 
and different impurity potentials when two impurities 
are nearest neighbors compared to impurities infinitely 
far separated from each other. The model presented be- 
fore will go astray if one places a large patch of impurity 
atoms next to a large patch of host atoms. In this case 
there will be a charge donation of the host into the impu- 
rity, independent of the distance between them, thus even 
if the impurity atom is very far away from the host. In re- 
ality this will cost Coulomb energy, which scales with the 
distance and should prevent too much of a charge flow. 
It is thus important to do self consistent calculations. 
Density Functional Theory captures these effects reason- 
ably well. The method, language and general features 
described in the previous sections are directly transfer- 
able to density-functional super cell calculations. 



A. Cluster shape dependence for Lio.75Mgo.25 

Mean-field approximated model calculations including 
several thousands of sites are possible and finite size ef- 
fects therefore can be largely eliminated. Density func- 
tional theory calculations are more involved and although 
modern computers are very powerful it is important to 
try to optimize cluster sizes and shapes. Naturally Mean- 
field theories like most implementations of Density func- 
tional theory are still numerically easier to solve then 
true many-body calculations, from which some knowl- 
edge about optimal cluster sizes can be obtained. In 
Fig. ([7]) we show three calculations for the LDA band- 
structure of Li4 (middle panels) and LiaMgi (top panels), 
calculated with the use of LMTO.^^ The calculations all 
correspond to different ordered structures of the same 
chemical composition. The difference between the three 
rows is the shape of the super cell. In the middle row 
one can see in orange that the folded band-structure de- 
pends on the choice of the super cell used. The unfolded 
bands are independent of this choice (when no impurities 
are introduced). In the top row one can see the unfolded 
bands when one Li atom is replaced by Mg. If one com- 
pares the bands of LisMgi with the folded bands of Li4 
then a striking feature appears. Looking for example to 
the lowest s derived band at -4 to eV going from F 
to in the unfolded Brillouin zone. Whenever the Li4 
folded bands cross the Li4 main bands (see middle rows) 
a gaping of the bands at that place in k space appears 
when one Li atom is replaced by Mg. 

This is strongly related to the difference between the 
calculation of an ordered structure and the calculation 
of a randomly substituted material. For ordered super- 
structures gaps appear when the folded bands cross the 
original bands. The dependence of this folding on the 
shape of the super cell opens the possibility to not only 
average over different configurations, i.e. different impu- 
rity densities and different orderings of the impurities, 
but also to average over different cluster shapes. 

The beauty of ah initio calculations over model calcu- 
lations is that one for each calculation assuming a certain 
impurity configuration and super cell shape, not only ob- 
tains the electronic properties, but also a total energy. 
This total energy can be used as a weighting factor when 
the different cluster sizes are averaged. On the bottom of 
Fig. ([7j) we show the total energy of the 3 different super 
cell calculations per site, with respect to the energy of the 
lowest energy state. The material studied will be grown 
or annealed at a certain temperature, where one can as- 
sume the system to be in thermal equilibrium. Then one 
could rapidly cool the system. The distribution of the 
different clusters can be assumed to be governed by the 
Bolzman statistics one would have at the annealing tem- 
perature. Annealing the LiMg crystal at 300K would give 
a distribution of 0.29, 0.33 and 0.38 for the three differ- 
ent super cells as shown in Fig. ([7]). If the sample would 
be in thermal equilibrium at 30 K the distribution func- 
tion would be 0.05, 0.23 and 0.72. Naturally the before 
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FIG. 8: (color online) LDA band-structure of Li75Mgo.25 cal- 
culated using a super cell of size 2x2x2 {Nt = 16) for panel 
(a) and (b) or a super cell of size 4x4x4 [Nt — 128) for 
panel (c) and (d). The calculations are averaged over differ- 
ent configurations, whereby in panel (a) and (c) the number 
of Mg atoms is fixed at 4 and 32 and only the random posi- 
tions varies between the different configurations, whereas in 
panel (b) and (d) the number of Mg atoms in the super cell 
[Ns) as well as their random positions varies. The intensities 
are scaled according to probabilities given by the distribution 
function in Eq.{|22]) with Nt = 16 {Nt = 128) and p = 0.25. 



sketched numbers are just exemplary. A real calculation 
needs much larger cluster sizes and the average over more 
then three different configurations or super cell shapes in 
order to give realistic results. 



B. Cluster size and averaging over the number of 
impurities in different super cells for Lio.75Mgo.25 



FIG. 9: (color online) LDA band- structure of Li with H sub- 
stitutions for different concentrations. The band-structure is 
calculated for a cluster of 128 atoms with random impurity 
configurations, unfolded to the original Brillouin zone. 



of the same Fig. ([8]). Here the cluster size is increased 
from 2x2x2 to 4x4x4. Again in panel (c) the total 
number of substituted sites is kept fixed at Ns = 32 for 
each configuration which varies the random positions of 
these substituted sites. In panel (d) also the number of 
substituted sites is varied between the different configu- 
rations. For larger clusters (at these impurity levels) it 
becomes less important to have the exact random distri- 
bution as given in Eq. (|22]) . One furthermore observes 
that for a 4 X 4 X 4 cluster of a bcc structure containing 
A^^ = 128 sites the averaged band-structure already looks 
very reasonable. 



The model calculations as shown in the previous sec- 
tion show averages over configurations where both the 
possible number of substituted atoms as well as the ran- 
dom placement of these atoms is varied. One would ex- 
pect that for large enough cluster sizes the variation of 
the number of impurities becomes less important. In Fig. 
dHI) we show in panel (a) a calculation done for several 
configurations where for each configuration the number 
of substituted atoms is chosen to be 4, which for a super 
cell of 16 sites results in a density of p = 0.25. In panel 
(b) we show a calculation for the same super cell, but 
now we also vary the total number of substituted atoms. 
The calculation in panel (b) includes averaging over all 
configurations according to their distribution function as 
given in Eq. (|22|) with Nt = 16 and p = 0.25. The 
average density is p = 0.25 in both cases. One can see 
that including an average over different number of substi- 
tuted atoms for small cluster sizes improves the calcula- 
tion. One can compare these results to panel (c) and (d) 



C. Concentration dependence in H substituted Li 

To illustrate the effect of disorder on the band struc- 
ture as obtained from density functional theory, we per- 
form DFT calculation of Li metal substituted with Hy- 
drogen using TB-LMTO code and super cell approxima- 
tion. The size of the super cell is 4 x 4 x 4 (128 sites 
due to bcc structure). A random number generator is 
used to generate 100 configurations for any given H con- 
centration. The unfolding procedure is used to present 
the calculated band structures in the original BZ. All 100 
band structures are plotted in the same plot. Fig. (j9j) 
shows such plots for 6.25, 12.5, 25 and 50% of Hydro- 
gen. The impurity effect is twofold, first, the host band 
structure is broadened in both energy and momentum, 
and, second, there is an upward shift in energy with re- 
spect to the Fermi energy. The shifts can be understood 
in terms of on-site energy differences between Li and H 
states. Since the ionization energy of H is about 8.2eV 
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(a) 





FIG. 10: (color online) LDA band-structure and density of 
states for Lio.5Ho.5- On the left calculations are shown for an 
ordered Li H structure. On the right calculations are shown 
for a disordered system. Calculations are done for a 128 site 
super cell with random impurity distributions, unfolded to 
the original Brillouin zone. 



higher than that of Li, the H 5-states are expected to 
have higher binding energy than that of Li. This means 
that the Hydrogen impurity is an acceptor. The shifts, 
however, are compensated somewhat by the changes in 
the band widths. At about 50 percent, the band struc- 
ture of LiH alloy turns over to that of pure H which is 
also broadened and shows the shifts in the Fermi energy 
which is now due to Li+ impurities as can be seen in the 
right bottom panel of Fig. ([9]). 



D. Order versus disorder in the Lio.5 Ho. 5 alloy 

The specific case of Lii/2Hi/2 is rather interesting since 
it allows us to study the effect of order on the electronic 
structure. The electronic structure of the ordered alloy 
is shown in the left columns of Fig. (p!Q|) . It is calculated 
with the conventional bcc unit cell that corresponds to 
two interpenetrating simple cubic lattices. The material 
is an insulator with the gap of about leV which separates 
the occupied H Is band with the conduction band formed 
by the states of mainly Li character. The random alloy 
is, however, metallic (see right columns of Fig. (p!Q|) ). We 
note the very large smearing effect of the band structure 



especially in the region dominated by lithium p based 
states above the Fermi energy. There is, of course, no 
background nor smearing of the bands in the ordered 
case. 

IV. CONCLUSION 

In conclusion we have discussed a band structure based 
procedure which can be successfully used to describe the 
electronic structure of disordered alloys and randomly 
substituted systems. We have used a simple tight bind- 
ing based model to describe the physics of the differences 
between minimum size super cell calculations and large 
super cells with random distributions of impurity atoms 
in the cell. In order to visualize the changes we present 
a simple unfolding procedure which is demonstrated to 
work very well indeed. We have shown that with small 
impurity potentials the band structure is qualitatively 
described by a rigid shift as in a virtual crystal approxi- 
mation. However even here the bands become broadened 
in both energy and momentum space. This broadening 
becomes very large for larger impurity potentials. We 
also show that the Fermi surfaces are strongly smeared 
out for medium size impurity potentials and that the 
peak position of the broadened structure as a function of 
momentum would yield a Fermi surface determined by 
the crossing at an energy corresponding to the electron 
count i.e. the "Fermi energy" which does not agree with 
the Luttenger theorem size of the Fermi surface. 

In order to quantify the effects of the impurity scat- 
tering effects we present the real and imaginary parts of 
the self energy which clearly demonstrates the various 
effects impurity scattering and randomness has on the 
electronic structure. The smearing effects at the Fermi 
energy can cause a semiconducting density of states in the 
ordered case turn into a metallic density of states. This 
effect is quite spectacular in the comparison of a density 
functional based ab initio calculation of the electronic 
structure of ordered and disordered LiH. The results we 
show for Mg Li and Li H alloys using density functional 
methods demonstrate that ah initio methods can be used 
which opens the door to the study of a large number of 
real material systems. 
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